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Abstract 



In this paper, a technique for constructing quasiperiodic structures 
is suggested, which allows one by the assigned matching to restore the 
atoms density distribution formula of a corresponding quasicrystal. 
The algorithm to restore the atom density distribution has been con- 
sidered on the example of the Penrose matching. The analytical record 
of a Penrose quasicrystal is givenQ. 

1 Introduction 

In the basis of traditional crystallography is the conception of periodically 
repeated structures. I.e. in crystals there is realized such a packing of one 
structural cell, which possesses translational symmetry in respect to a shift by 
some vector. Translational symmetry is a definition of a long-range order in 
crystals. Simultaneously with translational symmetry, rotational symmetry 
through angles 



can also exist. The same situation was in the studies on tilings, packings, 
mosaics etc.: it was believed that if there existed not one, but two or several 
matching elements, one could find a method of periodic tiling of a plane with 
the elements. Such tilings also possess rotational symmetries (|1|). However, 
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in 1974 R. Penrose [1] showed that with several elements, two as minimum, 
one could match a plane quasiperiodically. Such quasiperiodic matchings 
can have symmetry axes, different from symmetry exes allowed for periodic 
structures. (The concept of symmetry for quasiperiodic structures is still 
to be defined). Penrose's mathematic discovery anticipated the existence of 
quasiperiodic structures obtained experimentally in metal alloys. In 1984 
an article was published [2], in which an experimental evidence was pre- 
sented for the existence of an alloy with exceptional properties. An X-ray 
photograph showed a sharp system of spots testifying for the existence of 
a pentagonal order (to be more exact, decagonal), and the availability of a 
long-range order. At present, a great number of papers have been published 
devoted both to experimental and theoretical studies of such structures called 
quasicrystals. There exist several methods to describe quasicrystals: a tech- 
nique based on projecting on a 2D or 3D space of a periodic hyperlattice 
of a higher dimensionality [3-6], another method of projecting [7] consists 
in constructing periodic packings of icosahedrons or pentagons in curved 3D 
and 2D spaces respectively, where symmetry exclusions are weaker than in 
zero curviture spaces. But a return to the real Euclide world requires to 
introduce disclination lines which generate a long-range icosahedral order for 
a 3D space or a long-rang order corresponding to the Penrose matching for 
a 2D space. In the works [8-11] a detailed analysis of the Penrose problem 
has been undertaken, the GDM (generalized dual method) rules formulated 
[8,9] or multigrid ones [10,11], allowing one to construct quasicrystals of an 
infinite extent with icosahedral symmetry (3D) and point symmetry of the 
fifth order (2D). A cluster approach to consider the quasicrystal problem has 
been suggested in [12]. In the work, a technique to describe quasicrystal 
structures is suggested, using a traditional approach to describe ideal (peri- 
odic) crystals and incommensurate structures. In particular the symmetries 
of quasicrystal structures are considered with the superspace groups [13,14], 
which now are a classical method of describing incommensurate structures. 



2 ID quasicrystals 

Even though there is no analogue of orientational order in one dimension, one- 
dimensional quasicrystals will be useful to illustrate some important prop- 
erties of quasiperiodic structures. To begin with, remember the principle of 
constructing periodic crystals. Consider a one- dimensional periodic lattice 
with the position of sites at the points 

Xn = n. (2) 
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In this case we have a periodic sequence of one interval of the length S = 1, 
as shown in Fig.[l|. To each point of the lattice we relate some group of 



Figure 1: The periodic sequence of one interval of the length S = 1. 

atoms — the basis, with all the groups being identical by the composition and 
arrangenment. The basis, consisting of J atoms, is defined by a J set of Xj 
coordinates which determine the positions of the atom centres with respect to 
a lattice point, to which the basis is related. Atoms constituting the basis are 
arranged with respect to a given lattice point in such a way that ^ Xj < S. 
A set of atoms, thus assigned, is a unit cell of a periodic structure. For such 
a structure, the atoms density distribution (we will consider atoms as points) 
p{x) can be written in the form of the sum 



p{x) 



X 



n 



Xj) , 



(3) 



where the first summation (j = 1,..., J) is taken with respect to all basis 
atoms, and the second - to all lattice sites. Quasiperiodic structures, in 
contrast to periodic structures, consist of several structural cells forming a 
quasiperiodic sequence. Structural cells of a quasicrystal can differ both by 
the number of atoms and their arrangement and by sizes. Consider a qua- 
sicrystal consisting of two structural cells of the same size, but with different 
bases. Denote the internal atom coordinates of the first and second bases 
of 

•^ii' -^^2 respectively. In order to obtain the atoms density distribution 
formula for a quasicrystal, analogous to (j^), introduce the function 



p{n + (3) 



n + l + (3 




n + p 


a 




a 



(4) 



where [xJ is an greatest integer part of the number, a > 1 is an irrational 
number, /3 is the arbitrary number. The function p(n+/3) (n = 0, ±1, ±2 . . . ) 
assigns a quasiperiodic sequence of zeros and units [10,15-17]. In the sequence 
the relation of the number of units to the number of zeros equals l/(cr — 1). 
Define the functions 



Pi(n + /3) 



l-p{n + (3) 
p{n + /3) 



(5) 
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as probabilities of arranging the first and the second bases in a cell at the 
assigned n. For such a structure the atoms density distribution p[x,(3) can 
be written in the form of the sum 



Jk 



p{.x,(3)= ^ ^'Pk{n + (3)^5{x -n- 



X 



3k, 



(6) 



-ao k=l 



jk=^ 



where ^ Xj^^ < S, the first summation {jk = I, .., Jk) is made in all atoms 
of the basis, the second summation {k = 1,2)- in all bases, and the third- 
in all sites of the lattice. In the case of a quasiperiodic structure with two 
intervals L and S between the neighboring sites, position of sites at the Xn 
points is assigned by the formula 



Xn{a, (3) = n + a + 



X 



n + p 



a 



(7) 



where x > 0, a is an arbitrary number, further we will consider (3 /a < 1, 
since in the opposite case a changes, i.e. a' = a + l/x[P/o'\- The rela- 
tionship describes the position of the lattice sites so that the interval 
Axn = Xn+1 — Xn between the sites has the property: 
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if pi{n + (3) 



L = 1 + - , if P2{n + (3) = 1 . 
X 



The parameter x determines the relationship between the interval lengths in 
the sequence. For instance, for a particular case of(T = x = T = (l + "\/5)/2 
the quasiperiodic sequence of two intervals is shown in Fig.|3. In this sequence 
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Fi gure 2; The quasiperiodic sequence of two intervals L and S. 

the interval L is repeated once or two times in succession, the interval S only 
once. For a quasiperiodic structures with two structural cells of different 
length, the atoms density distribution p{x,a,/3) can be written in the form 
of the sum 

oo 2 i/fc 

p{x,a,(3)= ^ ^pk{n + (3)^6{x-Xn{a,(3)-XjJ, (9) 

n=— oo fc=l ifc=l 
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where ^ Xj^ < S, ^ Xj^ < L. To further study structural properties 
of quasicrystals, make a detailed analysis of the functions Pkin + (3). To 
demonstrate, consider the function of a continuous variable x 



fix, 61,62 



X + 61 




X + ^2 


a 




a 



(10) 



where o ^ 6\> 62- The function /(x, 6x^62) is periodic in x with the period 
a. The graph j{x, 6x^62) represents a sequence of impulses with the length 
of {6\ — 62) and the distance between the impulses (u — 5i + ^2) (Fig.|]). For 




X 



2a 



3a 



4a 



Figure 3: The graph /(cc, (^2) represents a sequence of impulses with the length of 
{b\ — 62) and the distance between the impulses (a — Si + 62)- 

f{x, 61, 62) the following equation holds true: 

f{x + e, 61, 62) = fix, 61 + e,62 + e). 

Moving along the of the graph of the function f{x,6i,62) to n from 

the initial position of x = (3 will give the value of the function 



f{n + P,6,,62) 





n + /3 + 6i 




n + (3 + 62 




a 




a 


Pi 


{n + (3), f{n 


+ (3,1,0) =p2{n 



'11] 



For example: f{n + P,a,l] 
the position of the point A(n + /5) within the period a when shifting from 
the origin of coordinates by the value of n + (3. To do this, calculate first by 
what maximum integer M of a periods the movement occurred: 



M 



n + P 



a 



Then (Fig.|) 



A{n + p) =n + p- aM = a 



n + f3 



a 



(12) 



(13) 
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o 2(5 3a Mo {M+l)o 



A(n+|3) 



p l+(3 2+|3 3 + p 4 + |3 5 + p n + p n+l+ p 

Figure 4: A(n + /3) = n + /? - crM = cr {(n + /3)/cr}. 

where {x} is a fractional part of the number, {x} = x — \_x\. Taking account 
of the periodicity A(n + the function f{n + P, 61,62) can be defined as 
follows: 



f{n + P, 61,62) 



0, if ^ A{n + p) < a ~ 61, a - 62 ^ A(n + p) < a, 

1, a a-6i^A{n + P)<a-62. 

(14) 



Such a definition of f{n+P, 61, 62), equivalent to the definition of this function 
in ([TT|), will be useful in studying 2D and 3D quasicrystals. Divide the period 
a into several sections of the length 

Si = a - 61, S2 = 61 - 62, . . . , Sk-1 = 6k~2 - 5k-i, sk = 6k~i, (15) 

where cr > 5i > ^2 > . . . > 6k-i- Now the number of structural cells of 
the same size, but with different bases, can be increased until K, defining K 
probabilities as follows: 

Pi(n + p) = f{n + p, a, 61), p2{n + p) = f{n + p, 61, 62), ... , 

PK-i{n + P) = f{n + P,6k-2,Sk-i), (16) 

PK{n + P) = f{n + P,6K-i,0). 

For such a structure, the atoms density distribution p{x,P) can be written 
in the form of the sum 

ca K Jk 

p{x,p)= J2 (1'^) 

n=-oo k=l ifc=l 

where ^ Xj^ < S (/c = 1,2, ... ,K). Repetition frequency of the basis 
with k index in such a structure will be equal to Sk/cr. In the case of a 
quasiperiodic structure with two intervals defined in formula (|), divide the 
period a into sections as follows: 



si = a - 61, S2 = 61- 62 
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where a > 6i > 62 > ■ ■ ■ > Sm~i > 1- 



^m+l 1 ^m+li • • • j (19) 

where 1 > 5m> > ■ ■ ■ > Sk-i > 0. Then the first m of the probabihties 
Pi{n+P), ..pm{n+P), defined in (pJ^, will correspond to the arrangement of m 
bases within the interval with the length S, and {K — m) of the probabilities 
Pm+i{n + /?),... ,pK{n + P) — to the arrangement of {K — m) bases within 
the interval with the length L. For such a quasiperiodic structure, the atoms 
density distribution p{x,f3) can be written in the form of the sum 

oo K Ji^ 

p{x,p)= ^ ^pfc(n + /?)25(x-x„(«,/3)-x,J, (20) 

n=-oo fc=l jk=l 

where ^ x^,, < S (/c = 1, 2, . . . , m), ^ Xj^ < L (fc = m + 1, m + 2, . . . ,K). 
To construct a quasiperiodic structure of an assigned configuration with bases 
alternating in a certain manner, take advantage of a famous fact from the 
probability theory. If pi,p2,... , Pm are the probabilities of setting in m 
independent developments, then 

m 

P = l[pk (21) 
fc=i 

is the probability of setting in all m developments simultaneously. To illus- 
trate the above, consider a particular example: x = c"" = and pi,p2,Xn 
are taken from formulae (^), (|^) respectively. The sequence of intervals for 
Xn{(y,l3) will correspond to Fig.^]. For the given quasiperiodic sequence of 
the intervals, introduce some structural cells with the lengths L and S with 
different bases (Fig.|]). As seen from Fig.^ structural cells with the length S 
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Figure 5: The quasiperiodic sequence of six structural cells with the lengths L and S. 

are arranged as follows: on the right of the cell Ai, in one interval, is a cell 
with the length S; on the left of the cell A2, in one interval, is a cell with 
the length S; on both sides of the cell ^3, in one interval, is a cell with the 
length L. Structural cells with the length L in a sequence are arranged as 
follows: Bi - on both sides there are two cells with the length S; B2, B3 - 
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are behind each other: B2 is the left, the right structural cells. According 
to the above said, define probabilities corresponding to these structural cells 



PAi{n + (3) 


= Pi 


[n + 


[n 


+ 2 + /5), 


(22a) 


PA2{n + f3) 


= Pi 


[n + (3)pi 


[n 


-2 + /3), 


(22b) 


PA3{n + f3) 


= Pi 


[n + (3)p2 


[n 


-2 + /?)p2(n + 2 + /?), 


(22c) 


PBi{n + /3) 


= P2 


[n + (3)pi 


[n 




(22d) 


PB2{n + /3) 


= P2 


[n + (3)p2 


[n 




(22e) 


PB3{n + /3) 


= P2 


[n + P)P2 


[n 


+ 


(22f) 



Calculations for equations (^) (refer to Appendix ^ produce the following 
result. 



PA2 in + /3 
PA,{n + (3 
PBi {n + (3 
PB,{n + f3 
PB3{n + /3 



n + 13 + T 



T 



n + l + l/r2 + /3 



n + /3 + 2/r 




n + /3 + 1 


T 




T 



n + p +1 + 1/1 



n + p + l/r 



n + p + 1 



T 



n + p + 2/i 



n + p + l/r^ 



n + p + 1/ 

T 



T 



n + (3+l/T 




T 






n + 13 






T 





(23a) 
(23b) 
(23c) 
(23d) 
(23e) 
(23f) 



These probabilities assign quasiperiodic sequences of zeros and units (n = 
0,1,... ,61;/3 = 0) 

PAi- 10000100000001000010000000100000001000010000000100001000000010 
PA2- 00000001000000010000100000001000000010000100000001000010000000 
PAi- 00100000001000000000000100000001000000000000100000000000010000 
PBi- 00000010000000100001000000010000000100001000000010000100000001 
PB2- 01001000010010000100001001000010010000100001001000010000100100 
PB3. 00010000100100001000010010000100100001000010010000100001001000 
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3 2D and 3D quasicrystals 



The orientational symmetry is a new element in extending a quasiperiodic 
translational order from one dimension to two or three ones. 2D and 3D 
quasicrystals allow symmetries excluded for periodic crystals. A simple gen- 
eralization of one-dimensional quasicrystals will be 3D quasicrystals with 
vectors of the main translations ao, ai, a2 of the following form: 

oo Ko , Ki , K2 

no, ni, 712=— 00 fco,fci,fc2=l 

J. . . ^24) 
where rn{a, f3) = aoXno + aix„, + ^2X^2, 

^jkok,k2 = ^^jkok,k2 + ai?/ifcofcifc2 + ^2%fc„fc,fc, is the vector determining the 
location of the basis atom centres with the index /co^i^2 relative to the lat- 
tice point to which the basis is related, = rii + ai + (l/xi) [('^^ + 
Atoms forming the basis are arranged relative to a given lattice site in such 
a way that ^ x,^^^^,,^ < So, ^ y,^^^^^^ < Si, ^ z.-^,^^^^^ < S2 
{ki = 1,2,... ,mi, 1 = 0,1,2, So = Si = S2 = l), 
^ ^Jkok,k2 < Lo, ^ y,^^^^,^^ < Li, < z,,^^^^,^^ < L2 
{ki = mi + l,... , Ki, / = 0, 1, 2, = 1 + 1/xi). 
As an example, consider 2D and 3D quasicrystals whose corresponding prob- 
abilities pi{n + j3), P2in + j3) are defined by formula (H). These quasicrystals 
will consist of four structural cells in the case of a 2D quasicrystal, and eight 
ones in the case of a 3D quasicrystal. Make a more detailed analysis of 2D 
quasicrystals. Introduce the probabilities 

Pa{n + (3) = pi{no + Po)pi{ni + Pi) , (25a) 

Pbin + /3) = piino + /3o)P2K + A) , (25b) 

Pe(n + (3)= PiijiQ + /3o) Piijii + /3i) , (25c) 

Pdip. + /3) = P2{nQ + /3o) P2(ni + A) (25d) 

and structural cells A, B,C ,D corresponding to these probabilities. Consider 
the probabilities as functions of continuous variables x, y. These functions 
are periodic in x with the period uo, in y with the period ui. To demonstrate, 
plot the graphs of these functions in the plane and point out the areas in 
which probabilities Pa,Pb,Pc,Pd assume values equal to unity (Fig.^. As seen 
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Figure 6: The areas Ua, Uj,, Uc, U^, in which probabiHties pa, pb, pc, assume values 
equal to unity respectively. 



from the graph: 

Pa(x,y) = l, if 

P6(x,?/) = 1, if 

Pc(a;,y) = l, if 



rioo-o ^ X < (no + 1) (To - 1, 
niOi ^ y < (ni + l)o-i - 1, 

nodo ^ X < (no + 1) CTo - 1, 

(ni + l)a-i - 1 ^ y < (ni + 1) (Ji, 

(no + 1) ctq - 1 ^ X < (no + 1) (To, 
^ y < (ni + 1) (Ti - 1, 



1, if 



(26a) 
(26b) 
(26c) 
(26d) 



(no + 1) To - 1 ^ a; < (no + 1) do, 
(ni + l)(Ti - 1 ^ y < (ni + l)ai. 

Taking account of the periodicity p^, p^, instead of the conditions (p6|) 

on the variables x, ?/ throughout the plane XY ^ one can consider conditions 
on the values 



Ao(x) 



Ai(y) 



X 



X 



X 



. y 

ai \ — 



(27) 



in the rectangle Oq x ox placed in initial of coordinate. 
Then formula (^6]) can be written in the form: 

Pa(x,?/) 

Pd{x,y) 



1, if 0^ Ao(x) <ao-l, 0^ Ai(y) <ai-l, 

1, if ^ Ao(x) < (To - 1, di - 1 ^ Ai{y) < , 

1, if (To - 1 ^ Ao(x) < (To, ^ Ai{y) < (Ti - 1 , 

1, if (To - 1 ^ Ao(x) < (To, (Ti - 1 ^ Ai(y) < (Ti . 



(28a) 
(28b) 
(28c) 
(28d) 
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The probabilites Pa{n + ^^(n + Pc(n + Pd(n + /3) and conditions 
for them, analogous to (pTf ), we will obtain at x = tiq + /5o, y = rii + The 
quasicrystal, thus assigned for a particular case of (Tq = o"i = r, is shown in 
Fig.^ In increasing the number of structural cells, analogously to (p!6D, each 
region in the graph Fig.^ will be divided into several rectangulars, within 
which, the corresponding probabilities will assume the value equal to unity 
(Fig.||). As well as in a one-dimensional case (formulae |2^, reciprocal 
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Figure 7: The quasicrystal, assigned Figure 8: In increasing the number 

by formulae for a particular case of of structural cells, analogously to (|l6|), 

(To = o"! = T. each region in the graph Fig.|^ will be 

divided into several rectangulars. 

location of the bases can be correlated. Using the properties of probabilities, 
described in Section ^, one can construct a sufficiently broad class of the 
assigned quasiperiodic structures, although in this case one may have to take 
account of an infinite number of products of probabilities, analogous to (p2|), 
correlating reciprocal location of the bases. In the above considered technique 
of describing quasicrystals for a definite cell, a corresponding probability had 
the form of a product of one-dimensional probabilities. It is not always that 
with such a presentation of probabilities one can describe quasiperiodic struc- 
tures. In a more general case, probabilities will be corresponded by regions 
of an arbitrary form restrained by curves (Fig.^). Analytical presentation of 
such quasiperiodic structures will be the following: 



Pfc(n + /3) 



1, if (AoK + /?o),Ai(ni+A)) G f/fc, , . 

0, if (Ao(no + /5o),Ai(ni+/?i))^t/fc, ^ ' 
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Figure 9: In a general case, probabilities will be corresponded by regions of an arbitrary 
form restrained by curves. 



oo K Jfc 

p(n,/3)= + (29b) 

no,ni=— oofe=l ife=l 

Such a technique of describing allows one by a quasiperiodic ornament with 
the assigned matching rule to restore all the probabilities of + To 
understand better, how by the assigned matching to construct a quasicrystal, 
consider a particular example. 



3.1 Example: the Penrose Matching 

To construct a quasiperiodic ornament consisting of two rhombi, whose sharp 
angles are equal to 36° and 72°, the rhombi are decorated with sections of 



straight lines [9] (Fig . [T0|) and a rule is assigned that in making mosaic, sec- 
tions on the rhombi are to form continuous lines. With such a matching 
rule a quasiperiodic Penrose ornament is produced, which intersects a penta- 
grid [9] (Fig.|IUp known as the Ammann quasilattice. Otherwise one can say 
that a Penrose ornament, constructed from rhombi, decorates the Ammann 
lattice or is a matching on the Ammann lattice. The position of parallel lines 
in the Ammann lattice is given by the formula [9]: 



Xni =ni + ai + - 

T 



ni + Pi 



T 



(30) 



Construct a Penrose matching on a beforehand assigned quasiperiodic lat- 
tice (i.e. a with the assigned Oq, cti, Po, Pi), consisting of a set of lines 
parallel to the vectors eg, ei so that the lines being formed from the sec- 
tions, decorating rhombi, should coincide with the assigned lines. In the 
Ammann lattice, intersecting the Penrose tiling, eliminate the lines parallel 
to 62, 63, 64. In Fig.|TT| such a quasiperiodic lattice, decorated with a Penrose 
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Co 



en= (cos {2%n/5), (sin i2%n/5)) 



Figure 10: Ammann line decoration of a portion of a Penrose tiling. The tiles have been 
decorated as shown at top to reveal the Ammann quasilattice associated with the tiling. 



mosaic, is shown. Each lattice cell can be seen as a mosaic element covered 
with a certain pattern. A quasiperiodic structure, corresponding to the Pen- 
rose matching, will be produced, if in each site of the matching one atom 
is located (Fig.|T^. Such a quasicrystal consists of 23 bases (Fig.^^. Note 
that 23 bases consist of as few as 12 atoms. One can find probabilities of 
each basis in two ways: The first way is to correlate reciprocal location of 
quasicrystals bases using the properties of probabilities multiplication, con- 
sidered in Section 1. The second way follows from the description principle 
itself of a quasicrystal. To define a probability of any basis, with the as- 
signed no and rii, calculate A(no + Po) and A(ni + Pi), then the set of the 
points (A(no + Po), A(ni -|- Pi)) for each basis will thickly fill a certain region 
in the probability distribution graph. In calculating by the both ways the 



probability distribution graph will be of the form shown in Fig.lT3 
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Figure 11: Ammann lattice, intersecting the Penrose tiling, with reserved lines parallel 
to eo, ei. In such a quasiperiodic lattice cell can be seen as a mosaic element covered with 
a certain pattern. 



Analytic presentation of the probabilities is of the form0 (|29|a). Thus, 
the atoms density distribution p(n, (3) for the Penrose quasicrystal can be 
written in the form of the sum 

oo 23 Jk 

no,ni=-oo k=l jfc=l 

Since in the Penrose crystal the locations of some atoms in different bases co- 
incide, one can find probabilities corresponding to each atom (altogether 12). 



In Fig.|T5| the regions are shown, in which the probabilties, corresponding 
to definite atoms, are equal to unity. It is clear that, if several atoms belong 
to one basis, then the regions of such atoms can overlap. In this inter- 
pretation analytic presentation of the atoms density distribution p(n, /3) can 
be written in the form of the sum : 

oo 12 

no, ni=— oo j=l 



^Complete analytic presentation of the probabilities is given in the Appendix 
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Figure 12: 23 bases of a Penrose quasicrystal consist of as few as 12 atoms. 
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Figure 14: The probability distribution graphs have been obtained with a numerical 
analysis of the Penrose quasicrystal. (a) 30 by 30 cells, (6) 50 by 50 cells. Points of one 
colour correspond to a definite basis. 
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Figure 15: The probability distribution graph of the Penrose quasicrystal, in which the 
areas Ufc, corresponding to definite atoms. 



Thus, on the axample of the Penrose mosaic the rules of constructing a 
quasicrystal corresponding to the given matching have been obtained. One 
can construct a quasicrystal by the assigned quasiperiodic matching only if 
the following conditions are available: 

(a) If there are several mosaic elements decorated with sections of straight 
lines, so that in constructing a quasiperiodic matching from mosaic elements, 
the sections of straight lines form a grid consisting of L groups of parallel 
continuous lines. 

(b) Position of lines in a general case is assigned by the formula: 

XnXo^uf^i) =ni + ai + — 

Xi 

where 1 ^ / ^ L, or 

=ni + ai (34) 

as this takes place, in the case of ( ^41) the ai parameters must be known. 
Then the order of constructing a quasicrystal by the assigned matching will 
be the following: 

(1) Reserve two groups of parallel lines from (^) or (P^), and assign the 



ni + Pi 



(33) 
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parameters ao, ai, /^o? A and also the main vectors of the translations ao, ai. 

(2) On the assigned straight lines construct a quasiperiodic mathcing, so that 
the lines being formed from sections, decorating mosaic elements, should ex- 
actly coincide with the assigned straight lines. 

(3) In the produced picture, eliminate the lines parallel to a2, a^, ... a^. 

(4) Correlate a quasicrystal atom with every site of the produced matching. 
Then we get a set of several bases. 

(5) In the produced (only graphically so far) quasicrystal, calculate a set 
of Ai{ni + Pi) (/ = 0,1) values for each basis. Subject to the condition 
that — A^"; ^ Ti; ^ A'^;, and A*"; ^ 1, a set of points with the coordinates 
(Ao(no + Po), Ai(?7,i + Pi)) for each basis will thickly fill a certain region in 
the probability distribution graph. 

To demonstrate, we have considered a constructing of a quasicrystal on the 
plane. Since computer modelling and analysis of three-dimensional lattices 
is rather difficult, restrict ourselves only to the description of a rule for con- 
structing three dimensional quasicrystals. In the three-dimensional space, 
three dimensional figures — mosaic elements. For example, two rhombohe- 
dra will be Penrose rombus analogues (Fig. |16]) [9]. 




Co 



Figure 16: Two types of rhombohcdra that arise in icosahedral quasicrystal packings. 
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Three-dimensional mosaic elements are decorated with sections-areas of faces. 
A generalization of a Penrose rule for a three-dimensional quasiperiodic 
matching will be the following: in constructing a matching, the areas of 
the planes, decorating three-dimensional mosaic elements, form continuous 
planes. The planes are to constitute several groups of parallel intersecting 
planes. Position of parallel planes is defined by the law ( p3D or (Q). In the 
three-dimensional case intersecting lines of the planes will be crystallographic 
exes. For such a quasiperiodic matching, a corresponding quasicrystal is con- 
structed according to the following rule: 

(1) Reserve three groups of parallel planes, whose position is determined by 
the law ( p3|) or (pi]), and assign values of the parameters ai, Pi and also the 
main vectors of the translations a; (/ = 0, 1, 2). 

(2) At given set of the planes, construct a quasiperiodic mathcing, so that 
the planes perpendicular to the vectors ag x ai, ai x a2, a2 x ag formed 
from the plane areas, decorating three-dimensional mosaic elements, should 
exactly coincide with the assigned planes. 

(3) In the produced volume picture eliminate the planes non-perpendicular 
to the vectors ao x ai, ai x a2, a2 x ao . 

(4) Correlate a quasicrystal atom with every site of the produced matching. 
In this case we obtain a set of several bases. 

(5) In the produced (only graphically so far) quasicrystal, calculate a set of 
values Ai{ni + Pi) (/ = 0, 1, 2) . Subject to the condition —Ni ^ rii ^ Ni, 
and A^"; ^ 1, a set of points with the coordinates (Ao(no + Po), Ai(r;,i + 
Pi), A2(n2 -l- P2)) for each basis will thickly fill a certain area of the proba- 
bility distribution graph. 



4 Difraction in quasicrystals 

Radiation intensity, scattered on a crystal, is described by the formula 

2 



/(q, a, /3) 



dr exp(— z q-r)p(r, a, j3) 



V 



|p(q, a, (3)\ 



(35) 



where V is the crystal volume, p(q, ck, (3) is the fourier-image of density dis- 
tribution of the atoms, q = ag go + a]; gi + sl^ q2, a* = a^ x afc/[aj-(aj x a^)]. 
To calculate a fourie-image of the atoms density distribution for a quasicrys- 
tal, present Xm defined by formula (^) as 



X. 



ni(l + 



1 



1 

Xi 



ni + Pi 
o-i 



A 



+ ai. 



(36) 
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Then the fourie-image of the atoms density distribution will be of the form 
p(q, a, /3)) = - Pk{+f3)llF{l,ni)Y,eMi^-rjJ, (37) 

na ni, n2 =-oo /=0 jk 

where F(Z, ni) = exp{i qinai + i qirji - iqi/xi{{ni + Pi)/(Ji}), 

-fi = 1 + 1/ (aixi), rii = [3i/{aixi) + ai. 
Consider in the right-hand side of (^) the following expression 

vl/(n + A q) = riexp (-z * ^xp (zq-r,,) • (38) 

Since the ^'(x, q) function is periodic in terms of x with the period ((Tq, o"i, (72), 
then \l/(x, q) can be expanded into a Fourie series in terms of x in the par- 
allelepiped with the sides do, ai, (T2. 

00 / 2 \ 

^(x, q)= c(m,q)exp -z27r^^ . (39) 

mo, mi, m2 =— oo \ /=0 / 

The Fourie coefficients c(m, q) are calculated as follows: 

c(m,q) = -^/ c/x^(x, q)exp I i27r) ^^^^ ) , (40) 



a- JVc 



[, q) exp U 27r 2^ — — I 



where is the parallelepiped volume with the sides (Tq, cti, (72. Calculation 
of the integral in the right-hand side of (|40|) is reduced to the calculation of 
K integrals from the regions Vk of the parallelepiped V^, in which the cor- 
responding probabilities Pfc(x) take the values equal to unity. Then formula 
(|40|) will assume the following form 

Jk 

c(m,q) = Y Cfc(m,q)exp («q-r^ J , (41) 

where 

Cfc(m, q) = / c/x exp I ix; > M — ] ) . (42) 



The fourie-image of the atoms distribution function will be the following 



p(q,a,/3) =^X]^^P 

n, m \ i=0 ^ 



27rm; 



X ^ Cfc(m,q)exp(zq-r^.J . 
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Summing in n obtain 



2 

p(q, a, /3) 5(q - q^^) — exp (i qir]i - i 2tx^^\ 

n,m 1=0 V '^i / 



K Jk 

Cfe(m, q) exp(? q-r^- 

fc=i A-=l 



(44) 



where 



E/ 2?™; 2Txmi 



i=0 



Using formula write the Fourie series of p{r,ct,f3) 

2 



p(r,Q;,/3) =^ JJexp ( 

n m n ^ 



n, m ;=0 

K Jk 

X 

k=l jk=l 



iqiTTji -i2tx 



^ ^ Cfc(m, q) exp [i q- (r + r^.J j 



(45) 



Then radiation intensity scattered on a quasicrystal can be written as 
/(q,Q;,/3) = :^ /" exp(zq-r)p(r,Q;,/3)p(q,Q;,/3) 

^2 



^ K,K'Jk,Jk' 

7^ XI XI XI Cfc(m,q)cfc/(m',q)5(q- q„m) 



n , m 



X 



exp [ i q- (rj^ - r^^J + z ^ 2n^{mi> - nii) 



X / c/rexp (2r-(qn'm' - qnm)) • (46) 
7y 



Since ^ ^ 1 the integral in the right-hand side of (^q) has the following 
values 



/ dr exp {i r 

■ (,qn'm' - qnm)) 

Jv 



V , if Qn'm' — qnm — , 
, if qn'm' - qnm 7^ . 



I.e. 
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(a) 



KqAO) 





I(q,0,0.87) 




Figure 17; Fouric-spcctrum of a quasipcriodic sequence of 6 functions, assigned by for- 
mula (0), (a) n = 7, (6) n = 17. 



/ exp r- (qn^m' - qnm)) = ^ n 

•^^ 1=0 



(48) 



Thus, the intensity of the scattered radiation is not equal to zero, when 
n = n', m = m'. 

I{q,(X,f3) = 

^ K,K'Jk,Jk' 

= V^Yl 5^ Cfc(™'q)cfc'(m,q)5(q-qnm)exp(^zq-(rj^-rj^,)j . 



n,m k,k' 



(49) 



The following result is produced: intensity of a scattered radiation on an 
infinite crystal is independent of the a and /3 parameters 



/(q,a,/3) = J(q,0,0) = J(q). 



(50) 
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Figure 18: The Fourie-spectrum (a) of a Penrose quasicrystal (6). 



5 Symmetry Transformations for 
Quasicrustals 

Since physical property of a quasicrystal, radiation scattering, is independent 
of ex. and f3, introduce an invariance condition for a quasicrystal: if under the 
operator influence of the transformation A on the atoms density distribution 
p(r, a, (3) only the ol and /3 parameters change, but not the p function itself 

Ap(r,a,/3) = p(r,a',/3') (51) 

then the quasicrystal is invariant in respect to such a transformation. 

5.1 Translational Invariance 

As has been shown above, a quasicrystal can be modelled with a quasiperi- 
odic or periodic lattice. Each site of the lattice is related to some group of 
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atoms — the basis, there being several such bases, and their repetition fre- 
quency is determined by the Pfc(n + /3) probabihties. Let us see whether the 
quasicrystal, thus modelled, is translationally invariant. Shift the lattice in 
such a way that the lattice m site should be in the point with the coordinates 
a, then in the lattice, produced as result of the shift, the coordinate of any 
point will be determined by the vector 

rn(a, /3') = r^+mict, (3) - r^(0, (3). (52) 

For each component we obtain the following result 



{ai, PI) = {ni + mi + ai + 



Xi 



ni + mi + /3i 



-{mi + 



Xi 



mi + Pi 



(Ji 



nl + al^ 

Xi 



ni_ ^ mi + Pi 

(Jl CTl 



m + Pi 

(Tl 



nz + H 

Xi 



ni 

— + 



mi + Pi 
o-i 



I.e. 



Xni{ai, Pi) = ni + ai + 



Xi 



ni + ai {{mi + Pi)/cri} 



■ (53) 



(54) 



The probabilities are the functions of the A;(n; + Pi) values for which the 
equation 

Aiim + mi+ pi) = Ai{ni + ai {{mi + A)M}) (55) 
holds true. Thus, in translating a quasicrystal get the following formula 

T(m)p(r,a,/3) = p(r,a,/3'), (56) 

where T(m) is the translation operator, PI = ai{{Pi + m;)/cr;}. Since the 
translation transformation for a quasicrystal satisfies the condition (pi]), a 
quasicrystal is invariant in respect to translation transformations. Consider 
an operation of translation multiplication 

T(m(2))T(m«)p(r,a,/3) = T{ui^^^)p{r,a, (3^'^) = p{r,oi, /3^'^^^^). (57) 

Write out the f3 parameter transformation for an individual component 
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Figure 19: Superposition of two Penrose mosaics: Si = 0.12/t, 82 = 0.18/r. 

, (1) _ J A + "^P^ 



For the {x} function the equation {x} = {a; + n} = ±1, ±2, . . . ) holds 
true, therefore the equation of the form {{x} + y} can be transformed as 
follows 

{{x] + y] = {{x] + \x\ + y] = [x + y] (59) 

it means that 

/?/^^^^)=a.| ^^ + "^'^| + "^''' |. (60) 

Then for translation operators we get the following condition 

T(m(2)) T(m(^)) = T(m(2) + m^^)) . (61) 
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Translation operators, corresponding to the property of (^) form a cyclic 
group [17]. If in shifting the quasicrystal 1^5; = |A/a; — + mi)/ai\, 
large areas of the initial and the shifted quasicrystals will coincide (Fig.^). 

5.2 Rotational Symmetries 

As has been mentioned earlier, quasicrystals can possess rotational symme- 
tries not allowed for ordinary crystals. Consider rotational symmetries on 
the example of the Penrose matching. Reserve in the Ammann quasilattice, 
intersecting the Penrose tiling, the lines parallel to the eo, ei, 64 vectors 
(Fig.[20|). The Penrose quasicrystal can be considered in the coordinate axes 




Figure 20: The Ammann quasilattice, intersecting the Penrose tiUng, with reserved hnes 
parallel to Gq, ei 64 vectors. 

xo, x\ (Fig.|l^), or xo, (Fig.^). In the both cases a crystal is produced 
with the same structural cells and the common law of their distribution, only 
the a;, (3i parameters will change. I.e. densities of the atoms distribution in 
the both cases have the form (PTj). Assign the values of parameters in the 
first case, since any lattice site can be placed at the origin of the coordinates: 
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Figure 21: The Ammann quasilatticc, intersecting the Penrose tihng, with reserved hnes 
parallel to eg, 64 vectors. 

ao = tti = 0, ^ Po, Pi < T. To find parameters for a quasicrystal considered 
in the exes parallel to eo, 64, select one of the cells, partially superimposed 
with the crystal cell, constructed in the axes, parallel to eo, ei and being at 
the origin of coordinates. After the cell has been selected, the parameters 
Oq, 04 are found geometrically. The parameters /Jq, j3l are calculated by the 
probability distribution graph (Fig.|T^). The first cell corresponding to the 
probability pk, will give a value region of the parameters /3q, PI within the 
region Uk'- On moving by a single step in the direction of Gq, or from 
the region Uk one can get to several Uk' regions permissible for realization 
of the Penrose quasicrystal. If on moving by a single step from the Uk re- 
gion the moving is possible only to a definite Uk' region, corresponding to 
a cell in the assigned quasicrystal, then the u^' region of value parameters 
will constrict, it is clear that u^' C Uk'- Correlating a cell, selected at the 
origin of coordinates with severel cells of the quasicrystal will reduce the u^ 
region of value parameters Pq, PI- The described technique to determine 
parameters can be presented as follows: consider a probability distribution 
graph of a Penrose quasicrystal, periodically continued throughout the whole 



27 



plane, and a two-dimensional square lattice with a single step. A Penrose 
quasicrystal in the axes, parallel to eg, 64 will be produced, if a lattice is 
superposed on the probability distribution graph in such a way that a point 
of the lattice with the coordinates (rio, "^.4) should be within the 11^' region of 
a corresponding structural cell (Fig.^2|) under the number (rio, ^4). Introduce 



y 




(3i T 2t 3t 4x 5x 



Figure 22: A Penrose quasicrystal will be produced, if a two-dimensional square lattice 
with a single step is superposed on the probability distribution graph, periodically con- 
tinued throughout the whole plane, than a point of the lattice with the coordinates (a;, y) 
should be within the Ufe region of a corresponding structural cell. 

a rotation operator i?((y9,ro), where is the rotation angle, rg is the point 
coordinate, around which rotation is performed , then for the Penrose crystal 

i?(27r/5, ro)p(r, a, /3) = p(r, a', (62) 

The considered quasicrystal is invariant in respect to rotations through an- 
gles, multiple of 72°, i.e. it has the symmetry axis of the fifth order. In a 
general quasicrystal will have the symmetry axis of the q order, if the 

following condition holds true 

i?(27r/g,ro)p(r,a,/3) =p(r,a',/3') (63) 

In such a quasicrystal the position of lines is assigned by formula (|30[) with 
Xo = Xi = • • • = Xfc, c^o = 0"! = • • • = CTfc, = a(cos(27rn/g), sin(27rn/g)). 
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The transformations (^) are a cyclic group of Cq rotations. A complete 
list of groups of proper rotations of a two-dimensional quasicrystal will be 
supplemented with a group of the same rotations taken together with 
reflections in respect to q axes. Three dimensional quasicrystals can also be 
invariant in respect to the T, W, P transformation groups, retaining invariant 
a regular tetrahedron, a cube (or octahedron) and a dodecahedron [9] (or an 
icosahedron) , respectively. 

6 Conclusion. 

To describe quasicrystals, in contrast to periodic crystals, use is made of 
a concept on qyasicrystals as objects consisting of several structural units 
having a quasiperiodic ordering. It is due to a quasiperiodic ordering of 
structural units that symmetries, different from classical, are permissible for 
quasicrystals. The developed approach to the description of quasicrystals 
allows simple analytic expressions to be obtained for the main structural 
characteristics — probabilities, which determine the repetition frequency of 
corresponding structural cells. A periodic crystal can be presented as a par- 
ticular case of a quasicrystal with one structural cell and a corresponding 
probability with the period 1. Such a technique of describing quasicrystals 
generalizes a definition of the crystal, but does not complete it, since there 
exists a mathematical model of aperiodic crystals^. It is not excluded that 
there may exist still more complex crystal structures. 

A Function product calculation 

of the form /(x, Si, 62) • f{x — e, 6^, 64) 

Function products of the form f{x, 61,62) ■ f{x — 6,63,64) are calculated 
graphically. For generality of results consider the f{x,6k,6k+i) functions 
whose impulse length is greater than the distances between the impulses 
6k - 4+1 > (T - 6k + 6k+i- 
(a) In the case (Fig.pSj), when e < a — 61 + 6^. 

f{x, 61, 62) fix - e, 63, 64) = fix, 63 - e, 62) (64) 

when (Fig.^ a > e > a - 61 + 62: 

fix, 61, 62) fix - e, 63, 6i) = fix, 61,64 + a - e) + fix, 63 - e, 62) (65) 

^to be published 
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f(x,bi,b2) 





f(x-e, 53,84) 





/(x,5i,52)/(x-8, 83,54) = /(x,53-e,52) 







Figure 23: e < a - Si + 5^ 



f(x,bi,b2) 





f(x-e, 53,84) 











f(x, 8i,82)/(x -8, 83, 84) = f(x,bi, 84 - o + £) + /(x, 83 -£,82) 







Figure 24: a>e>a-Si + 52 
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(b) li e > a, /(x, 53,54) is shifted more than by one period. Make the 
following transformation: 

f{x - e, 53, Si) = f{x - a{e/ a} - a [e/ crj , ^3, ^4) = 

/(x - (T{e/(T}, 53,54), (66) 

after that one can consider the case (a) and (b) for a shift by the value 
a{e/a}. 



B Analytic presentation of probabilities 
for a quasicrystal defined of formula ^ 

A 7^ (32- 



Pi(n + /3) = 1, if 



P2(n + /3) = 1, if 



P3(n + /3) = 1, if 



P4(n + /3) = 1, if 



P5(n + /3) = 1, if 



P6(n + /3) = 1, if 



P7(n + /3) = 1, if 



A(no + /5o 
A(ni 

A(no + /3o 
A(ni 

A(ni+/3i 
A(ni+/5i 

A(no + /5o 
A(ni+/5i 
A(ni+/5i 

A(ni+/5i 
A(ni+/5i 
A(ni 

A(ni+/5i 

A(ni 

A(ni 

A(no + /9o 
A(ni 

A(no + /3o 
A(ni 

A(ni+/3i 
A(ni+/?i 



^ 1/r, A(ni +/5i) < r, 

> A(no + /3o) + 1/r, 

^ 1/r, A(ni +/3i) < r, 

< A(no + /So) + 1/r, 

> A(no + /3o)/r + l/r, 
>rA(no + /5o)-l/r2, 

^ 1/r, 

< A(no + /5o)/r + l/r, 
>rA(no + /5o)-l/r2, 

< r, 

> A(no + /9o)/r-l/r, 

< rA(no + /3o) - l/r^ 

> A(no + /5o), 

< A(no + /5o)/r- 1/r, 

< rA(no + /3o) - l/r^, 

< r, A(ni +/3i) ^ 1/r, 

< A(no + /5o)-l/r, 

< r, A(ni +/5i) ^ 1/r, 

> A(no + /3o)-l/r, 

< A(no + /3o)/r + l/r3, 

>rA(no + /3o)-l, 



(67a) 



(67b) 



(67c) 



(67d) 



(67e) 



(67f) 



(67? 
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P8(n + /3) = 1, if 



P9(n + /3) = 1, if 



Pio(n + /3) = l, if 



Pii(n + /3) = 1, if 



Pi2(n + /3) = 1, if 



Pi3(n + /3) = 1, if 



Pi4(n + /3) = 1, if 



Pi5(n + /3) = 1, if 



Pi6(n + /3) = 1, if 



Pi7(n + /3) = 1, if 



A(ni 
A(ni 

AK + /3o 
A(ni+/3i 

A(ni 

A(no + /?o 

A(ni 

A(ni 

A(no + /5o 
A(ni 

A(no + /9o 
A(ni 

A(ni+/3i 

A(no + /3o 
A(ni+A 
A(ni+A 

A(no + /3o 
A(ni+A 
A(ni+/3i 

A(no + /3o 
A(ni+A 
A(ni 

A(no + /9o 
A(ni 

A(ni+A 
A(ni 



< A(no + /3o)/r + l/r3, 
>rA(no + /5o)-l, 

< ^, 

> AK + /3o)/r + l/r3, 
<rA(no + /5o)-l, 

> A(ni+/?i), 

> A(no + /3o)/r + l/r3, 

>rA(no + /3o)-l, 

^ 0, A(ni + < r 

> A(no + /3o)/r + 2/r, 

< A(no + /3o)/r + 2/r 

> AK + /3o) + l, 

^0, 

>rA(no + /?o) + l/^ 

< AK + /5o) + l, 

< ^, 

<rA(r2o + /3o) + l/^ 

> A(no + (3o) + 1/r, 

< A{no + f3o) + 1/r, 

> A(no + /3o)/r + l/r, 

^ 1/r, A(ni+/5i) ^ 1/r, 

< AK + /3o)/r + l/r, 

< 1/r, A(no + /3o) ^ r, 
>rA(r2o + /3o)-l, 
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Pi8(n + /3) 



Pi9(n + /3) 



P2o(n + /3) 



if 



if 



if 



P2i(n + /3) = 1, if 



P22(n + /3) = 1, if 



A(ni 
A(ni 

A(ni+/3i 
A(ni 

A(ni+A 

A(ni 

A(ni 

A(ni 
A(ni 
A(ni 

A(ni 



<l/r, 

> A(no + /3o) -1/r, 

< rA(no + /5o) - 1, 

< 

< AK + /3o)-l/r, 

< A(no + /3o)/r-l/r2, 

> A(no + /5o)-l, 

< A(no + /3o)/r-l/r2, 

< A(no + /9o) -1, 

> rA(rio + /5o) -2, 

^ 0, A(r2o + /3o) < r, 

< rA(no + /5o) -2, 



(67r) 



(67s) 



(67t) 



(67u) 



(67\ 



In other cases the probabilities PA;(n + /3) = 0. 
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